Carex pulicaris abundance is positively associated with soil acidity, rainfall and floristic diversity in the eastern distribution range

Carex pulicaris is considered an endangered species, and further losses are forecast under the influence of even moderate climate change. Local studies indicate that temporal declines in C. pulicaris abundance are positively correlated to decreases in precipitation and increases in air temperature. Determining ecological properties on larger scales than local ones can help develop effective protection programs for the species. We hypothesize that the local relationships observed between C. pulicaris abundance and precipitation, air temperature and soil properties will be confirmed in a spatially-oriented large-scale study performed in situ. Therefore, the present study takes a novel, large-scale integrated approach to (1) precisely characterize the ecological requirements of C. pulicaris within its eastern distribution range, and (2) determine the influence of its community type, soil properties and climatic conditions on its abundance. It was found that C. pulicaris is not a dominant or codominant species in the studied phytocoenoses in the eastern distribution range. Five natural vegetation groups including C. pulicaris, with significantly diverse species compositions, were resolved: well supported Estonian, Polish, Slovak and Radecz groups, and a weakly-supported Ambiguous group. The abundance of C. pulicaris was found to be positively correlated with the composition of the geographically-diversified plant communities and atmospheric precipitation, and to be also negatively associated with latitude and soil pH. Although the species is adapted to a relatively wide range of soil types, such adaptation requires appropriate substrate moisture level and light conditions. The species prefers moist organic and mineral soils and grows on both acid and neutral medium, characterized by a narrow C:N ratio, with various amounts of digestible total P, Mg and N, and low levels of digestible K. Climate change, manifested by reduced rainfall, may be one of the most important predictors negatively affecting the occurrence of C. pulicaris.

www.nature.com/scientificreports/ phytosociological data (17 relevés), other phytosociological relevés of Poland (13 relevés), which are recorded in the Polish Vegetation Database 32 , as well as previously published data from Slovakia (13 relevés) 29 . In total, our analyses were based on 43 phytosociological relevés recorded from 27 localities (1-3 per locality). The syntaxonomic approach to communities was adopted after Matuszkiewicz (2006) 25 . The habitat conditions were assessed based on Ellenberg indicator values (EIVs) 33 : ecological (L-light value, T-temperature value, K-continental value, F-humidity value, R-acidification value, N-nitrogen value, Ssalinity value), sociological (Gr-group of vegetation classes, K.1-class, o-order, v-alliance, u-suballiance) frequency and risk (m-frequency of appearance, D-domination, A-the trend of change, G-threats). In the case of 17 of the examined locations, the soil data were also analysed.
Climatic data. Climatic conditions were evaluated for all the analysed locations based on the rasters collected in the https:// world grids. org library, developed for the needs of the "Soilgrids" project, organized by the International Soil Reference and Information Centre (ISRIC, https:// www. isric. org/ explo re/ soilg rids). The data are currently available in the archive at https:// web. archi ve. org/ web/ 20170 61905 4443/ http:// www. world grids. org/ doku. php/ start.
The analyses included latitude, longitude, daily minimal temperatures, daily mean temperatures, daily mean temperatures in December/January, and annual mean monthly precipitation ( Table 1). The raster resolution for temperature data was 1 km, and 5.6 km for rainfall data.
Chemical composition of soil samples. Composite samples were collected from the surface layer of the rhizosphere (0-20 cm) for each of the natural sites 34 (Table 1). In the soil material, the following analyses were performed: loss on ignition (organic matter) by burning soil samples in a muffle furnace at the 550º C; total C, N and S content by elementary analysis (Costech Elementary Analyzer ECS 4010, Italy); pH in H 2 O and pH in 1 mol dm −3 KCl was determined potentiometrically; salinity by conductometry.
In addition, the content of Mg, K, Ca, Na, Cd, Co, Cu, Ni, Pb, Mn, Fe and Zn soluble in 0.5 mol dm −3 HCl (the so-called available forms) 35,36 and soluble in concentrated HNO 3 and HClO 4 , at the ratio 1:1 (the so-called total forms) was determined using an ICE series 3000 spectrometer with flame atomization (FAAS) 37 . The content of Na, Ca and K was determined by flame atomic emission spectrometry, while that of the other elements was determined using flame atomic absorption spectrometry. The limits of detection were (mg kg −1 ): Ca-0.004; Mg-0.002; K-0.001; Na-0.004; Fe-0.004; Cd-0.003; Co-0.010; Cu-0.005; Ni-0.008; Pb-0.013;  Table 1. Colouring of symbols: black-Estonian group, blue-Slovakian group, orange-Radecz group, purple-Ambiguous group, red-Polish group, white-locality in Slovakia with the occurence of C. pulicaris, but lacks species composition data. www.nature.com/scientificreports/ Zn-0.003 and Mn-0.002. The content of available and total P was determined with spectrophotometric molybdenum blue method (690 nm wave length) 37 using a Marcel MEDIA™ spectrophotometer. The accuracy and precision of the analytical methods and procedures used were confirmed using certified reference material: CRM036-050 Loamy Sand 4 (CRM 036-050 produced by Resource Technology Corporation, USA and UK). The effectiveness of the process was validated with 90-95% efficiency. The results shown are the mean values of three measurements, with working standards made from Merck standards at a concentration of 1000 mg dm −3 .
The soil type classification follows the international and regional standards 38-40 . Statistical analyses. The abundance of Carex pulicaris was compared with the phytosociological composition of sites and three sets of environmental variables: mean Ellenberg indicator values (EIVs), climate variables (CLIM) and soil variables (SOIL). The species composition and soil types of the sites were compared by multivariate analyses using the vegan package in R environment: non-metric multidimensional scaling (NMDS) of the species composition data matrix (43 sites × 310 species), and principal component analysis (PCA) of the soil data matrix (17 sites × 34 soil www.nature.com/scientificreports/ variables) 41,42 . The raw species composition data matrix and Bray dissimilarity index were subjected to NMDS, and raw data were centred and scaled prior to PCA. Cluster analysis (normalised data, Ward agglomerative method) was performed to delineate natural phytosociological groups; the grouping pattern was used throughout results. The relationship of C. pulicaris to site structure and its relationship with possible explanatory variables were assessed by fitting C. pulicaris abundance data over NMDS or PCA ordinations using the envfit function of the vegan package 41,42 .
In addition, species composition on the sites was compared with environmental data. The collected data were utilised to construct six data matrices; three types of environmental data matrices: (1) mean Ellenberg indicator values (EIVs), (2) climate variables, (3) soil variables, and three corresponding species-composition data matrices. The species composition data in each locality were averaged per locality to obtain a reduced dataset corresponding to database-mined climate data (reduction from 43 sites to 27 localities); this dataset was independent of the species composition in each location.
Canonical correspondence analysis (CCA) was performed in vegan 41,42 to formally test for the ecological interpretation (EIVs and CLIM) of the sites and the species similarity assemblages. ANOVA was applied for statistical testing, and the global CCA model was statistically significant. The parsimonious CCA model was selected by a permutation test using forward model selection.
SOIL data were not subjected to CCA, since the number of sites (n = 16, with available plant community and soil variable data) was much lower than the number of soil variables (n = 34). Instead, it was tested whether soil data may explain species composition data; briefly, the envfit function was used (in vegan 41,42 ) to superimpose SOIL variables over NMDS ordination of species composition data (n = 16 sites in total).
To characterise each of the localities, plant diversity indices, i.e. the number of species recorded per relevé, Shannon's H diversity index, and Pielou`s Evenness were also calculated in the vegan package of R or in Past 3.10 41-43 . One-way analysis of numerical variance (ANOVA) and Tukey's HSD pairwise multiple comparison were used to test for differences between means of particular parameters of different vegetation groups. Prior to analysis, assumptions of ANOVA, data normality and homogeneity of variance, were tested applying the Shapiro-Wilk test and Levene's test, respectively. Slight deviations from data normality were tolerated. The level of significance was set at p < 0.05 in all analyses. Statistical analysis of the obtained results was performed using R or in Past 3.10 42,43 .
Regarding the tested chemical properties of the soil, the statistical significance of the differences between means was determined by testing the normality of distribution in each group and homogeneity of variance in all groups, followed by ANOVA with Tukey's post hoc test. The significance was set at p < 0.05. These analyses were performed using Statistica 12.5 PL software (StatSoft Inc., Tulsa, OK, USA).

Results
Vegetation types with the occurrence of C. pulicaris. In total, 310 moss and vascular plant species were recorded (59 mosses, 251 vascular plants). Carex pulicaris locations were found to demonstrate a range of species compositions. Our cluster analysis of the similarity in species composition between sites ( Fig. 2) revealed five natural vegetation groups, as presented throughout the manuscript. These groups roughly correspond to the region of origin: Estonian sites (Estonian group), Polish sites including Danish sites (Polish group), Slovakian sites (Slovak group), relevés from Radecz in Poland, which appear as a separate group (Radecz group) and a final group composed of sites of different regions (Ambiguous group); not surprisingly, the final group has low statistical support (Fig. 2).
Mosses are prevalent in the Radecz (11-15 species per site, mean 12.3) and Slovak groups (7-18 species per site, mean 11.3); however, they are rare or even absent in the rest, where number of moss species per site ranged from 0 to 8.   Table 1. Colouring of clusters: black-Estonian group, blue-Slovakian group, orange-Radecz group, purple-Ambiguous group, red-Polish group. www.nature.com/scientificreports/ within the Slovak group, those of the Ambiguous group are scattered across ordination, while those of the Estonian group are found on the left and lower part of the plot. The EIV indicators suggest that the left-right separation of relevés derives mostly from the presence of species with higher values for Light, sociological class and Humidity (all positively associated with CCA1), as well as Continentality, sociological group of vegetation classes, trend of change and Nitrogen (negatively associated with CCA1). The relevés are separated across the second CCA2 axis by soil reaction and frequency of appearance. The abundance of C. pulicaris is also significantly correlated with the CCA-constrained ordination of plant communities. In general, NMDS and CCA correspond to each other and the abundance of C. pulicaris is significantly associated with vegetation type, whose differentiation seems to be driven by factors acting on large geographic scales.
The NMDS analyses of the species composition data of the reduced dataset (27 localities, Fig. 5A) found the overall groupings identified in the NMDS results of the whole dataset to be preserved. Latitude, Longitude, Mean precipitation and Minimal daily temperature, as well as C. pulicaris abundance, were found to significantly correlate with NMDS ordination.  www.nature.com/scientificreports/ The CCA analyses performed on the independently-gathered database of climate data described above (Fig. 5B) found that while the reduced CCA model was statistically significant (p = 0.001), only the first two axes were statistically significant (p < 0.05). The forward selection procedure suggests that only Latitude and Mean precipitation are statistically significant variables (p < 0.01). However, this CCA model explains only 11.24% (constrained) of the total variance. The resulting general grouping pattern of vegetation relevés appears slightly changed, with the Estonian group and Polish relevés being found in the right and upper parts of the plot, and the Slovakian group in the lower left part (Fig. 5B). The abundance of C. pulicaris was found to fit significantly to the groupings observed in the reduced CCA model; it also correlated negatively with Latitude, which is reflected on CCA ordination (Fig. 5B).
Soil properties at sites with C. pulicaris. The studies on soils associated with C. pulicaris in Estonia, Poland, Denmark and Slovakia indicate that the species prefers moist habitats: peat bogs with various degrees of decay (Murshic Histosols) and humus mineral soils (Umbric Gleysols).  www.nature.com/scientificreports/ In Poland and Estonia, some sites were located in fens with a weakly moorsh forming process. At the surface level (rhizosphere), they were characterized by low ash content (< 20%) and acidic or slightly acidic pH. In Slovakia, C. pulicaris also grew in fen, but strongly decayed bogs with high ash content (> 40%) and slightly acidic or neutral pH. The rhizospheres of the bogs differed significantly with regard to total nitrogen content ( Table 2). The soils from Estonia contained on average about 80% more total nitrogen than those in Poland and Slovakia. No significant differences were found in mean total sulfur content, C/N ratio or salinity. The soils on which C. pulicaris occurred were characterized by a varying content of available phosphorus (from small to medium amounts) and magnesium (from small to very large amounts) and a very low content of available potassium. The Estonian, Polish and Slovak bogs did not differ significantly in their available K or Na content. The peatlands from Poland contained significantly more available P than those from Slovakia, which were richer in available Ca and Mg (Table 2). www.nature.com/scientificreports/ Carex pulicaris was also found on moist mineral soils (Umbric Gleysols) with a mean organic matter content of 10 to 16% and a C/N ratio ranging from 13 to 14. Despite their similar systematic affiliation, the surface soil layers (rhizospheres) differed in their chemical composition between Estonia, Denmark and Poland ( Table 2). The soils in Estonia were slightly acidic, with a low content of assimilable phosphorus and potassium but with an abundance of Mg. On the other hand, soils from Denmark and Poland were acidic and poor in assimilable K and Mg; however, the soils from Estonia and Denmark demonstrated significantly lower assimilable P and considerably lower exchangeable Na content than those from Poland. Moreover, all three groups of soils were characterized by a low salt concentration and a similar content of exchangeable Ca. All analyzed soils contained small amounts of bioavailable and total forms of heavy metals (Cd, Co, Cu, Pb, Ni, Zn, Mn and Fe; Tables 3 and 4).
The PCA analysis of the soil dataset (17 sites × 34 soil variables) found that the site groupings reflect generally geographic locations and soil types; in addition, this patterning appears to be irrespective of plant community type, and of the abundance of C. pulicaris (Fig. 6A,B). In contrast, the results of the NMDS analysis of plant communities with regard to soil subset (16 sites × 194 species, plant community data are missing for single locality from Slovakia-SK_PodR) are roughly in line with the relevé grouping in the whole dataset (43 sites × 310 species). NMDS analysis separated the Estonian and Polish vegetation groups, on the right of the plot, from the Ambiguous group and one Slovakian relevé, on the left of the plot (Fig. 6C). The only soil variables significantly fitted to this NMDS grouping were pH and total Pb. Most importantly, C. pulicaris abundance fits significantly to the NMDS groupings (p < 0.01) (Fig. 6C). Table 3. Content of available heavy metals in soils. www.nature.com/scientificreports/

Discussion
The Estonian, Slovak and Radecz groups, which were most closely related to a particular region within the range of C. pulicaris, were characterized by a homogeneous type of vegetation. The Estonian group was represented by communities of the Molinion caeruleae alliance, sometimes with a significant share of Carex davallianae, and the other two by communities of the Caricetalia davallianae order. In Slovakia, C. pulicaris is associated with non-forest communities of minerotrophic fens, particularly with the Sphagno warnstorfiani-Tomenthypnion alliance 29 . In turn, the Polish and Ambiguous groups showed greater diversity: their composition included communities from both of the above-mentioned syntaxonomic units, with the Polish group including some from the Caricetalia nigrae order, with some tendency towards the Nardo-Callunetea class. Throughout its entire range, C. pulicaris is commonly derived from the communities of the Molinion alliance, and this can be seen in Slovenia, the Netherlands and Sweden 7,44,45 . In France, it was recorded in acidic fen grasslands (Caricetalia nigrae communities) 46 , and in Great Britain, in communities such as calcicolous grassland community Sesleria albicans-Galium sterneri in the Carex pulicaris-Carex panicea subcommunity, mire community Carex dioica-Pinguicula vulgaris and Molinia caerulea-Cirsium dissectum fen-meadow (according to Table 4. Total content of macroelements and heavy metals in the soil. www.nature.com/scientificreports/  Table 1. Colouring corresponds to clustering pattern, symbol size reflects Carex pulicaris abundance. Carex pulicaris abundance is not related to PCA grouping patterns (p > 0.05) but significant to plant community grouping (p < 0.05). Colouring of symbols: black-Estonian group, blue-Slovakian group, orange-Radecz group, purple-Ambiguous group, red-Polish group, white-locality in Slovakia with the occurence of C. pulicaris, but lacks species composition data. Note: Averages from three subsamples were inputs for localities SK_Pod1 and SK_PodR. www.nature.com/scientificreports/ the UK NVC 47 ). The last of these communities is the British equivalent of Cirsio-Molinietum Sissingh & de Vries 1942, described from the Netherlands, Belgium, Germany and Ireland 48 . The groups of plant communities separated as a result of the analyses were generally characterized by considerable plant species diversity, with the only exception being the Estonian group (Fig. 3); the lower diversity indices in this group was undoubtedly influenced by the high share of Molinia caerulea in the studied phytocoenoses. This is a highly competitive grass, and its presence in large clumps not only limited the number of safe places for the germination of other plant species diaspores, but also hindered their development. This was especially true for annual plants, and to a lesser extent, the perennials, which show vegetative reproduction. Moreover, most of the plant patches were partially shaded by the crowns of trees growing in the immediate vicinity. Under these conditions, C. pulicaris often grew in the Molinia clumps. We observed this phenomenon not only in the locations in Estonia, but also in the Rozwarowskie Marshes in Poland (PL_BagR).

The localizations of Carex pulicaris
Regardless of the research area, C. pulicaris was often accompanied by species associated with grassland communities, particularly meso-and eutrophic hay meadows and riverside herbaceous plant communities; these are permanently, or at least periodically, moist (order Molinietalia caeruleae). The significant share of these species in some of the examined patches of C. pulicaris communities from the Caricetalia davallianae order (e.g. in Slovakia) may indicate the presence of fluctuations in the water regime and/or a slightly increased availability of nutrients for plants. Such coexistence of species may take place over a longer period of time, and in the case of a permanent change in habitat conditions, it may lead to another stage of succession towards meadow communities from the Molinion caeruleae alliance 29 . Hence, it is sometimes difficult to unequivocally define the syntaxonomic affiliation of plant communities from C. pulicaris.
The share of C. pulicaris in the analyzed groups was diversified, and it was neither a dominant nor a codominant species in any of the phytocoenoses. In contrast to the eastern end of its range, C. pulicaris can be found abundantly in the remainder of its distribution area; it has been reported as the dominant or codominant species in some phytocoenoses, for instance in Great Britain and Ireland 47,49,50 . As C. pulicaris is a poorly-germinating 51 , light-requiring species that is weak in competition, its abundance is favored by extensive land use, which limits the development of tall, expansive perennials. Our findings indicate that the mowing of meadows in areas where C. pulicaris occurs (Domysłów and Rozwarowskie Marshes) favors its survival, even under conditions of high competition from Molinia caerulea, and in the vicinity of Phragmites australis plantations. It is classified as a species dependent on mowing in the area of the fenoscandian limestone forest meadows 52 , and was only found in mowed areas in the costal grasslands in the Stockholm archipelago 45 . Also in the Netherlands, this type of land use is necessary to maintain the structure and species composition of meadows belonging to the Molinion caeruleae alliance, where C. pulicaris is present 7 . This prevents the succession towards scrub and forest, especially when the habitat is fragmented.
Our findings also indicate that a higher abundance of C. pulicaris is correlated with greater precipitation, albeit to a small extent (Fig. 5). This relationship is demonstrated by the fact that the species was most abundant at the Polana Biały Potok (Tatry Mts) site (PL_PolB), which had the highest mean annual precipitation (915 mm) of the examined sites. Carex pulicaris is a sub-Atlantic ranging type 53 , indicating that it favors areas with a significant amount of rainfall. It is a common species in Great Britain, being sometimes codominant in plant communities, and is often observed in densely-tufted patches with Sesleria albicans 47 ; in this country, it has been found in areas with an annual rainfall of 1239 mm 54 .
An important factor limiting the range of C. pulicaris, as well as other sub-Atlantic species, may be the occurrence of winter frost 55 . Comparing the isotherms of the mean month temperatures of the coldest month (°C) calculated for sea level 55 with the species distribution data 13,53 , a clear relationship can be seen between the limit of its range and the -4 °C isotherm. The results of local studies from mountain areas in southern Germany 14 highlight the significant role played by both temperature and precipitation on the occurrence of C. pulicaris. In contrast, while our present findings do not show any clear relationship between species abundance and air temperature, the CCA analysis (Fig. 5B) suggests that rainfall is of greater importance.
Our results suggest that in the eastern edge of the C. pulicaris distribution range, its abundance may be latitude dependent, with a lower abundance observed in plant communities in more northerly locations. In addition, it grew sparingly in the north-eastern part of our research area, including almost all the Estonian sites, but had a fairly significant share in the community patches in the south-east, including most Slovak sites. However, no such relationship was observed in the north-western extremities of its range, where in some plant communities, it is very common and is the dominant species 49 ; it should be noted though that, being a sub-Atlantic species, the conditions for development in these areas are very favorable.
Carex pulicaris is found on various types of bogs, including fen, transitional and lagg of raised bogs, as well as wet heathlands, Nardus grasslands and wet meadows [56][57][58] . In addition to organic soils, it also often inhabits humus and moist mineral soils 44,56 , while it is less common in other habitats, such as clay or sandy peat soils 59 . In the sites in the present study, C. pulicaris grew in humid but diversified habitats: from bogs, where the Murshic Histosols soils developed, to humus mineral Umbric Gleysols soils 38 , regardless of whether it was at the edge or the center of the distribution range. Carex pulicaris occurs on a variety of soil types with organic matter content ranging from 1.6 to nearly 100% 49 .
Certain chemical and physical properties of soils have a strong influence on the occurrence of plants 60 . Our research showed that the surface layers (rhizospheres) of the soils in the tested areas demonstrated great variation in their physical and chemical compositions. Murshic Histosols were acidic to slightly acidic, while Umbric Gleysols were acidic to neutral 38 .
The areas across Europe where C. pulicaris occurs demonstrate similar variations in soil pH 44,49,57,58,61 . However, C. pulicaris generally preferred slightly acidic soils, and did not favour strongly acidic and alkaline soils. Our study showed that soil pH was significantly correlated with the grouping of plant communities with C. pulicaris participation. Soil acidity affects such features as the trophic nature of the habitat, the quantity and quality of www.nature.com/scientificreports/ humus, and the abundance of assimilable macro-and microelements. The wide range of soil pH values tolerated by C. pulicaris may indicate a well-developed adaptability to various environmental conditions; however, this contrasts with previous observations, i.e. that the species occurs in bogs of a narrow ecological amplitude 29 . Nitrogen, K and P levels are particularly important for plant community stability. In wet habitats, an increase of P, K and N levels in soil can lead to a decline in the diversity and abundance of plant species 44,[62][63][64] . The results of our analyses showed that both organic and mineral soils were generally low in available P and K, which may support the growth of C. pulicaris. These habitats have previously been found to have low nitrogen content 44,65 , and data from areas within the C. pulicaris range suggest that its communities require low N, P and K concentrations 44,61,62,[65][66][67][68][69] ; however, our present results do not indicate any significant relationship between the abundance of C. pulicaris itself and the P, N and K content of the soil. The plant communities of C. pulicaris are sensitive to excess N-NH 4 levels and N-NO 3 deficiency; this relationship is related to long-term stagnation of water on the soil surface, resulting in the development of anaerobic conditions 70 .
Our own studies indicate that C. pulicaris abundance is not significantly affected by soil Ca, Mg and salt content. The species was recorded in soils with various Mg, Ca and salinity levels 44,65 , and its associated plant communities have been observed in habitats relatively rich in nutrients 29,49,65,67 . It has also been found to grow in soils containing low levels of heavy metals 39,40 .
Carex pulicaris appears quite resistant to anthropogenic influences known to increase habitat trophy, as indicated by our finding that macro-and micronutrient content appears to have no clear effect on its abundance in plant communities. On the other hand, distinct changes in soil moisture and the cessation of pratotechnical practices may limit its abundance 44,70 .
Carex pulicaris responds well to surface floods in the autumn-spring period as well as to strong falls in groundwater level in summer 71 . Short periodic floods cause a reduction in the mineralization of organic matter; they also drive the release of nutrients, which are beneficial for the habitat, and stimulate the formation of Fe, Ca and P complexes, which are normally inaccessible to plants 72 . However, if maintained for too long, excessive soil moisture levels result in changes in the oxidoreductive and chemical properties of the soil that are unfavorable for the growth of C. pulicaris 70 . In contrast, a rise in the level of groundwater may favor paludification processes, which in turn will lead to the disappearance of C. pulicaris as a result of the reconstruction of plant communities with more competitive members. In addition, excessive dehydration and oxygenation of the habitat causes rapid mineralization of organic matter and the release of large quantities of mineral nutrients; this is also unfavorable for C. pulicaris, and favors the entry of expansive perennials, trees and shrubs.
Finally, most of the Estonian sites were characterized by large shares of trees in the vegetation patches; this can increase shading and thus may be the reason for the low abundance of C. pulicaris. The specific humidity requirements described above explain the sensitivity of this species to changes in humidity in the habitat.

Conclusions
Within the eastern end of its range, C. pulicaris is present in five natural phytocoenological groups: four closely correspond to geographic regions, viz. Estonia, Poland, Slovakia and Radecz (Poland), and a fifth Ambiguous group scattered throughout the study region. In general, the groups show a significant diversity of species composition. The Estonian group was represented by communities from the Molinion caeruleae alliance, the Slovak and Radecz groups by communities from the Caricetalia davallianae order, while the other two showed greater diversity with regard to their syntaxonomic units. However, C. pulicaris was not dominant or codominant in any of the groups with regard to the studied phytocoenoses.
Our findings indicate that the abundance of C. pulicaris is positively correlated with the composition of its geographically-diversified plant communities and atmospheric precipitation. In the eastern part of its range, C. pulicaris grows on organic (Murshic Histosols) and mineral (Umbric Gleysols) soils associated with high groundwater conditions. These soils are characterized by an acidic to neutral reaction, narrow C:N ratio, low salinity, mostly low levels of available P and K and micronutrients, and varying levels of available Mg, Ca and Na. Carex pulicaris appeared to grow most abundantly in plant communities which prefer slightly acid soils.

Data availability
All data matrices that support the findings of this study are available from the corresponding author upon request.